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Abstract 

Diirr and Hoelbling recently observed that the continuum and chiral limits do not commute in 
the two dimensional, one flavor, Schwinger model with staggered fermions. I point out that such 
lack of commutativity can also be seen in four-dimensional staggered chiral perturbation theory 
(SXPT) in quenched or partially quenched quantities constructed to be particularly sensitive to the 
chiral limit. Although the physics involved in the SXPT examples is quite different from that in the 
Schwinger model, neither singularity seems to be connected to the trick of taking the n*'^ root of 
the fermion determinant to remove unwanted degrees of freedom ("tastes"). Further, I argue that 
the singularities in SXPT are absent in most commonly-computed quantities in the unquenched 
(full) QCD case and do not imply any unexpected systematic errors in recent MILC calculations 
with staggered fermions. 

PACS numbers: 12.39.Fe, 12.38. Gc, ll.lO.Kk 
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I. INTRODUCTION 



The one-flavor Schwinger model in two dimensions (2D) has recently been studied by 

n 

Diirr and Hoelbling 1] using both overlap and staggered fermions. The square root of 
the staggered determinant is used to eliminate the extra taste degree of freedom. This is 
jasicallv the same trick that is employed in current MILC three-flavor QCD simulations 
^1 Q 0] (the fourth root is required in 4D). Although many of the results in Ref. 
are encouraging for the use of "rooted" staggered quarks, there is one disturbing feature: 
The chiral and continuum limits of do not commute. If the continuum limit is taken 
first, the known continuum result is reproduced. But if the chiral limit is first taken for the 
staggered quarks, the continuum limit then disagrees with the exact result. This leads to 
three key questions: 

• Do similar singularities appear in 4D QCD with staggered fermions? 



arities induce uncontrolled systematic errors in the 



• If they do appear, do such singu 
results reported in, e.g., Refs. j2| 

• Are the singularities the result of the rooting procedure? 

Here, I address these questions for 4D QCD in the context of staggered chiral perturbation 
theory (SXPT) RBBIU Q. In SXPT. quantities fo. wh.eh the chi.a. and 

continuum limits fail to commute can easily be found in the quenched or partially quenched 
cases. As explained in Sec. [HI the way the noncommutativity comes about is simple: Taste 
violations split the masses of mesons, and the splitting becomes the dominant effect in the 
meson masses as the chiral limit is approached. When the physics is particularly sensitive 
to the chiral regime, the chiral and continuum limits do not commute. Quenching or partial 
quenching leads to such sensitivity through double poles in neutral meson propagators, 
which enhance the infrared (IR) regime in loop diagrams. In the context of SXPT, the lack 
of commutativity has nothing directly to do with the issue of taking the root of the staggered 
determinant, and occurs in normal, "unrooted," staggered theories also. 

In Sec. mil I discuss the issue of commutativity of limits for unquenched SXPT ("full 
SXPT"). Defining safe quantities as those for which the chiral and continuum limits com- 
mute, I argue that most standard physical quantities are safe because the full theory has no 
double poles and is therefore better behaved in the IR. In this sense it is very similar to full 



continuum chiral perturbation theory (XPT). The safe quantities include those commonly 
determined in simulations, such as /tt or the ratio of squared Goldstone meson mass to quark 
mass. However, the shorthand statement "standard quantities are safe" is misleading, since 
is in fact simple to write down other quantities for which the limits do not commute in full 
SXPT. A rough guide is that quantities that are finite in the chiral limit of full, continuum 
XPT are likely to be safe in this sense, although this guide itself has exceptions. Further- 
more, there is at this point no proof, but only an intuitive argument for why the singularities 
are absent to all orders for the standard quantities in the full case. However, as discussed 
in Sec. IIV| even if one assumes that the argument breaks down and singularities actually 
appear at some higher order, the associated new errors in previously computed quantities 
would be negligible. 

Finally, in Sec. I mention two other examples of lack of commutativity, one with 
Wilson fermions and one in the quenched 2D Schwinger model. I point out that the known 
examples seem to indicate that the rooting procedure and the commutativity issues are 
independent. I conclude with a discussion of whether noncommutativity might be seen in 
topological and r( physics in 4D. 



II. EXAMPLES FROM SXPT 



The simplest example of noncommutativity of limits in SXPT occurs in the quenched pion 
mass, where the "pion" is made of two degenerate valence quarks of mass m^. In continuum 
quenched XPT, we have at next-to- leading order (NLO) |l2l |: 

f Ml \ 1 -2ml , . ,.^ 

where A is the chiral scale; m-Q, the contribution to the rj' mass from the anomaly; /, the 
decay constant normahzed so that / ~ 131MeV; and Xwi the tree-level pion mass squared, 
Xvv = ^firriv (yU is a constant with dimensions of mass.) The ■ ■ ■ represents less singular 
terms that do not lead to noncommutativity of limits. Here, and throughout this paper, I 
work in the infinite volume limit for simplicity. 

The corresponding expression for the Goldstone pion in SXPT is Ql 



f Ml \ 1 -2ml , / r ,.x 



(2) 
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where Xvv is the squared tree-level mass of the taste-singlet pion, 

xiv = Xvv + a?^i = 2fimv + o?^i (3) 

with a^A/ the splitting of the taste-singlet pion from the Goldstone pion. In Eq. Q the 
neglected terms include effects of the taste- violating hairpins ||7|. These are suppressed by 
a? and cause no problems with the limits. 

Taking the ratio of lattice to continuum results, we have, to this order, 

- 1 + 1 z5!!o ,„ f 1 + . . . . (4) 



This ratio clearly goes to 1 if we take the continuum limit, a — ^ 0, first. But it blows up if 
the chiral limit (m^ — > 0, Xw ~^ 0) is taken first, because the staggered theory gives the 
wrong answer in this case. Of course XPT breaks down once the correction term gets large. 
In this case the breakdown occurs because of the "quenched chiral log" in the continuum. 
The logarithm is cut off in the staggered theory by the taste splitting, a^Aj. 

This quenched example makes my basic point, but it is not very similar to the effect 
seen in Ref. There, the staggered theory is found to give a finite result for either order 
of the limits rriy and a — > 0; it is just that the result is incorrect (different from the 
continuum) when my — is taken first. An example that has this kind of behavior can be 
found in the partially quenched 4D theory.^ 

For simplicity, I take degenerate sea quarks of mass in the partially quenched 
theory, and again consider a pion made from degenerate valence quarks of mass rriv. In the 
continuum, the pion mass at NLO is ^ 

i^).,. = '^Ti^i - ^ ■ ■ ■ • 

where Xss = 2/ims, and ■ ■ ■ represents analytic terms. 

In the staggered case, consider Np degenerate dynamical staggered fields with mass mg. 
1 define Ng here as the number of sea quarks in the continuum limit: Ng = Np if the fourth 
root of the determinant is taken; while Ng = 4Np if it is not. The result has the same form 
when written in terms of Ns, whether or not the root is taken. Putting sea and valence 

^ But note that the example in Ref. occurs in the full 2D theory. 
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quarks separately degenerate in Eq. (48) of Ref. [3], one easily arrives at the one-loop mass 
of the Goldstone pion: 

{^)^^ -'^y^^l - J (•^'>-/^) + ■ ■ ■ • 

where x^ss defined analogously to Xvv^ Eq. Q, and ■ ■ ■ represents analytic terms and the 
effects of taste-violating hairpins, which again cause no problem with the limits. 

If the continuum limit is taken first, Eq. © reproduces Eq. (jSJ, so the two will give 
identical results no matter how the chiral limit is subsequently taken. On the other hand, if 
first we take the valence chiral limit (m^ 0, with fixed), the lack of commutativity of 
limits exactly parallels the quenched case, with Xss playing the role of mg. As before, 0{a?) 
terms in the staggered theory cut off a chiral log that is divergent in the continuum theory. 

For my purposes, a more interesting chiral limit occurs when and mg both approach 
0, but with nis vanishing much more slowly, so that 

Xss ~ 1 . (7) 
ln(xw/A) 

as Xvv 0, with C a positive constant. In this limit (with a fixed in the staggered case) 
Ml \ 12 



where ■■■ represents additional 0{a^\a{a?)) or O{o^) terms, coming from taste-violating 
hairpins or taste- violating analytic terms. If we now take the continuum limit, 

" Ml 



(10) 



in disagreement with Eq. (jH)). On the other hand, the staggered theory clearly reproduces 
Eq. (jHJ if the continuum limit is taken first in Eq. (jU)). 

Note that the noncommutativity of limits has nothing to do with the fourth-root pre- 
scription, and has exactly the same form with or without the fourth root. Of course, for 
fixed number Np of dynamical staggered fields, the effective number of continuum sea quarks 
Ns depends on whether the root is taken, so we could say in that sense that the rooting 
prescription trivially affects the the strength of the singularity, but not its existence. 



5 



III. COMMUTATIVITY IN FULL SXPT 



The lack of commutativity of limits appears in the computed one-loop SXPT quantities 
only in the quenched or partially quenched cases, which are enhanced in the IR 
relative to the full case and therefore especially sensitive to the chiral limit. This is due 
to the double pole structure that shows up in flavor-neutral meson propagators. In the full 
QCD case, there are only single poles, and 1 therefore do not expect noncommutativity to be 
generic in this case. For example, if we first go to full QCD {itiy = nig = m; Xw = Xss = x) 
in Eqs. (0) and ©, the NLO corrections then always vanish in the subsequent m — » 0, a — > 
limits, independent of their order. 

The absence of the commutativity problem in most full QCD quantities seems to be a 
general feature of SXPT. It is equivalent to the statement that the limit m ^ at fixed a is 
smooth in full SXPT, i.e., that there are no IR divergences on shell in this limit. I believe 
this to be the case because, taking the additional limit a —>■ 0, one recovers the massless 
continuum XPT. Ordinary chiral power counting in this limit shows that IR divergences 
are absent: The derivative couplings suppress IR contributions and make loop effects finite 
even in the presence of massless particles.^ If we now turn on the lattice-dependent ((9(a^) 
and higher) vertices of SXPT, these act as effective mass terms and should not induce IR 
divergences that were previously absent. Power counting applied to SXPT then implies that 
any logarithms of that appear in the m — limit are "protected" by powers of a^, so 
they cause no problem in the subsequent a — limit. 

Full SXPT in the chiral limit is in fact closely analogous to full continuum XPT with 
massless up and down quarks but a nonvanishing strange quark mass. In both cases there 
are terms giving mass to some, but not all, of the pseudoscalar mesons: The pions remain 
massless in the continuum example; while the taste-^5 (Goldstone) meson remains massless 
in SXPT. And in both cases the mass or mass-like terms give rise to interactions lacking 
enough derivatives to suppress all IR divergences automatically. Yet, one expects no IR 
divergences in the continuum the strange quark mass turns on, and therefore I 

expect SXPT to be likewise well-behaved. 

However, a detailed proof that the chiral limit of full SXPT is "safe" may be rather 

^ Note that one needs to go on shell so that derivatives acting on external lines also suppress the diagrams 
in the massless limit. 
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delicate. Although the lattice-dependent vertices give no nonderivative interactions purely 
among the taste-^5 mesons j6[[7,Q (the analogous statement is also true in the continuum 
example), potential problems may arise when external non-Goldstone mesons (taste different 
from ,^5) scatter to give massless .^5 mesons in internal loops. 

As a simple example, consider massless SXPT theory with two flavors, and take only the 



3ory_v 

m 



coefficient C4 in the taste- violating, O(a^), potential [6|,|7[ to be nonzero. This is sufficient to 
give mass to all the non-Goldstone mesons. Now consider the contribution shown in Fig.^to 
the scattering of two taste-singlet, flavor charged, mesons through the exchange of taste-^5 
mesons. If the two vertices come from the C4 term, then there are no derivatives on the 
internal lines to suppress the logarithmic IR divergence. On shell, the momentum transfer 
through the loop will, in the generic case, provide a cutoff. However, with C4 vertices, there is 
an IR divergence in this diagram in the forward direction, pi = p^, p2 = p^. The divergence is 
canceled by the diagrams with one or both of the vertices replaced by kinetic energy vertices. 
The terms in the kinetic energy vertices where the derivatives act on the external lines 
provide the cancellation, possible because pi = ps implies pi ■ ps = —a^Aj = —QAa^C^^/p. 
Terms where the derivatives act on the internal lines are clearly not IR divergent; neither 
is the corresponding "s-channel" diagram (on shell). Though the result of this example is 
positive (no IR divergence on shell), it also shows some of the issues that would need to be 
addressed in a complete proof that full SXPT is safe. 

Pi P3 




FIG. 1: A "t-cliannel" scattering diagram for the scattering of two taste-singlet mesons, with 
internal taste-^5 mesons: T^f{pi) + ttJ {P2) —>■ T^fiPs) + t^J (Pi)- 

A further subtlety is the following: When I suggest that full SXPT is "safe," I really 
only mean that certain standard quantities are safe, i.e., they have commuting chiral and 
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continuum limits. A sufficient number of derivatives (with respect to quark mass) of any safe 
quantity will certainly produce IR divergences and hence lack of commutativity of limits. 
Which quantities are likely to be safe? The above arguments suggest that quantities that 
are free of IR divergences in the chiral limit of full continuum XPT will probably be safe in 



p 



m 



T. This means, for example, that pseudoscalar masses and decay constants as computed 



will be safe. On the other hand, a quantity like the pion charge radius,'^ which is defined 
in terms of an external current with nonvanishing momentum transfer, is IR divergent in 
the continuum chiral limit |14:]. The charge radius is therefore almost certain to have the 
wrong chiral limit at fixed a in SXPT; it should indeed behave much the same way as the 
quenched M^/2fimy, Eqs. (0) and 0. 

Furthermore, since in the continuum one can differentiate the squared meson mass once 
with respect to the quark mass and still have a finite quantity in the chiral limit, I expect 
that the corresponding quantities in SXPT are safe. For the squared Goldstone (taste ^5) 
mass this is fairly obvious from the exact lattice Ua{^)', one can also divide by (instead of 
differentiating with respect to) the quark mass, as in the full SXPT version of Eq. (jHl). It 
is less obvious that one can differentiate the squared masses of non-Goldstone (taste other 
than ^5) mesons and still have a safe quantity. For example, how do we know that a term like 
a^xlnx cannot appear in the squared non-Goldstone masses? Such terms are not divergent 
as m ^ 0, but their derivatives are divergent, which would lead to noncommutativity. A 
one-loop calculation, similar to the scattering calculation described above, shows that such 
terms do in fact cancel for the taste-singlet mass. Whether this really persists beyond one 
loop remains to be seen. 

Note that dividing a squared non-Goldstone mass by quark mass, rather than differen- 
tiating, is certainly unsafe, as is already clear at tree-level, Eq. Q. This shows why it is 
difficult to write down a simple rule that distinguishes safe from unsafe quantities in all 
cases. On the other hand, if we limit ourselves only to Goldstone mesons on the external 
lines, then it appears that SXPT quantities will be safe if the corresponding quantities in 
full continuum XPT are finite in the chiral limit. 



^ I thank S. Sharpe for this example 
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IV. ERRORS OF EXISTING NUMERICAL COMPUTATIONS 

Diirr and Hoelbling ^ worry that the possible noncommutativity of hmits with staggered 
fermions would induce new, uncontrolled, errors in the previously reported results with 
MILC staggered configurations 3,0,0,151. To the extent that SXPT describes the MILC 
simulations, we can argue that such errors are not present — even in the absence of a detailed 
proof that full SXPT is safe to all orders. That is because the extraction of physical results 
requires extrapolation only to the physical light quark masses, not to the chiral limit. The 
factor ln(x/A^) is of order 3 or 4 at the physical point {x ~ 140 MeV). Even should the 
logarithm not be "protected" by powers of x in some higher order of SXPT and represent 
a divergence in the chiral limit at fixed a, it would still be suppressed, numerically, by the 
additional powers of at that order. By looking at the numerics of the calculations in 
Re.. BBta, ,t ,s .ot Ka..d to co.v.ce oneself tKat a,, new en.. wouM ,nost .,.e>v be 
significantly smaller than the already-quoted systematic errors. 

As an example, consider the case of /^r. Imagine that there is in fact an "unsafe" contri- 
bution at NNLO (next-to-next-to-leading order), i.e., of order a^ln(x/A^). To estimate its 
size, define the dimensionless taste-violating chiral expansion parameter x^a by 

^ ' (11) 

where a^A is a "typical" taste- violating term at (9(a^). Taking for a^A the average pion 
splitting, gives ~ 0.09 on the MILC coarse lattices and ~ 0.03 on the fine lattices. 

Then I would expect the putative NNLO unsafe contribution to be generated by a term 
of size 

^J2^r,XB/A\ (12) 

B 

where the sum is over all 16 pion tastes, labeled by B. The unsafe contribution comes 
from the Goldstone pion in the sum. Note that I assume that the Goldstone pion appears 
only through such an average over tastes. That seems very likely, since the diagrams that 
distinguish separate tastes are diagrams with disconnected internal pion propagators, which 
at NLO occur just for the singlet, axial, and vector tastes. At NNLO there could be discon- 
nected internal Goldstone lines, but it is known that such contributions come with an 
extra factor of p^, making them safe. Putting in Inx/A^ ^ 4 and the values of x^a gives an 
unsafe contribution of ~ 0.2% on the coarse lattices and ~ 0.02% on the fine lattices. These 
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are clearly much smaller than other systematic errors (which are ~3% for /^r and ~1% for 
Ik/ fir) , even without taking into account the fact that some of the error would be removed 
by the continuum extrapolation (linear in Q;|a^ in Ref. [^). 

Of course, there would also be systematic errors coming from contributions of the other 
taste pions in Eq. (|12p. However, these contributions are safe, with commuting chiral and 
continuum limits. The associated errors are therefore standard higher order SXPT errors, 
and as such are already included in the error estimates of Ref. ^ . 

Note that my error estimate assumes that the "worst case" would be the existence of a 
singularity in full SXPT that would interfere with the extrapolation to the physical light 



quark masses 
in Refs. U 



his is reasonable despite the fact that intermediate stages in the analysis 
often involve fits to partially quenched lattice data, where we know that 
noncommutativity of limits does occur. The point is that the relevant NLO formulas such 
as Eq. (jH)) are known. We fit lattice data to them in a range of valence and sea masses that 
are far from the singularities. Since the light quark masses in the data are always more 
than a factor of three larger than the physical masses, the logarithms and the nonlinearities 
associated with them are always significantly smaller than at the physical point. It is 
therefore not surprising that (as checked directly in the calculations) the NLO corrections 
in the fits are under control and of the expected magnitude. If unexpected new singularities 
were to occur in higher order, the resulting errors in the partially quenched fits would be 
smaller than those induced in the subsequent full SXPT extrapolation, simply because the 
logarithms are smaller. 

Of course, if the staggered simulations had some additional noncommutativity of limits 
not captured by SXPT — having to do, say, with the fourth-root procedure — then an 
uncontrolled systematic error might in fact be present in the MILC results. However, recent 



advances in understanding the rootingprocedure 
Q to SXPT predictions flQ, 



30, In 



17|, coupled with the good fit of 



MILC data ^ to SXPT predictions [IlISi, make such a possibility seem increasingly unlikely. 

V. ADDITIONAL REMARKS 

The lack of commutativity of the chiral and continuum limits is not exclusively a property 
of staggered fermions. For Wilson fermions, it appears at one loop in a power counting for 
which 0{m) and O^a^) are treated as comparable At this order, a term of the form 
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a^X^T^X is present in the squared meson mass in the full QCD case.'* (The corresponding 
partially quenched case has not to my knowledge been studied.) After differentiating with 
respect to (or dividing by) the quark mass, such a term leads to an IR divergence as m — >■ 
for fixed a; whereas it clearly vanishes if the a limit is taken first. The source of problem 
is the absence of an exact chiral symmetry in the massless limit. This means that massless 
Wilson pions can interact without derivative couplings leading to IR singularities and 
failure of commutativity. 

Note that the Wilson case is significantly worse than the staggered case: In the latter, 
self- interactions of the massless Goldstone (taste ^5) pion are always proportional to at least 
two powers of momentum. At finite a, the non-Goldstone staggered pions are not required 
to have derivative interactions, but neither are they massless. 

My SXPT examples and certainly the Wilson case suggest that lack of commutativity 
has little or nothing to do with the rooting procedure. Is this also true of the 2D Schwinger 
models studied in Ref. The physics of the one-flavor 2D Schwinger model (ipip), where 
the noncommutativity is found, is quite distinct from the multifiavor 4D chiral theories 
discussed above. With one flavor, the "condensate" comes about because of the symmetry is 
violated by the anomaly, not because of spontaneous symmetry breaking. In the continuum, 
the anomaly leads to exact zero modes of the Dirac operator, which in turn saturate {ipip)- 
Since staggered fermions have only near zero modes in complex pairs, their effect cancels in 
the m — > limit, and {ipip) vanishes for fixed a. If the continuum limit is taken first, however, 
the complex pair of near zero modes becomes a pair of degenerate exact zero modes, and 
the rooting prescription works as desired, producing a single exact mode. (See Refs. jjj, |20| 
for details.) Since the physics here is inextricably tied to having only one flavor, however, 
it is not possible to separate cleanly the issues of the lack of commutativity and the rooting 
procedure within this model. 

In the two-flavor 2D Schwinger model, on the other hand, the rooting procedure is not 
required, yet I suspect that it would be possible to find noncommutativity in the masses 
of the "quasi-Goldstone bosons," of a similar nature to that seen in quenched or partially 
quenched 4D SXPT (Sec. |n|). That is because integration over a single (boson) pole in 2D 
has the same IR behavior as integration over a double pole in 4D. The problem here is 



I am grateful to O. Bar for pointing out to me this result in Ref. |18l |. 
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that there is no true condensate in 2D, so there is no simple chiral theory: {ipip) and the 
quasi- Goldstone boson squared masses vanish as fractional powers of m as m ^ 21]. It 
is therefore a nontrivial problem to discover an analogue of SXPT that would allow one to 
calculate the discretization effects caused by staggered fermions. Thus I am not able at this 
point to make a specific proposal for where to look for noncommutativity in this model. 

A simple alternative to the two-flavor theory in this context can be found in the quenched 
2D Schwinger model, which does not require the rooting procedure either (trivially). Ref- 
erence studies the quenched theory with staggered fermions, but since {^ip) blows up 
as 1/m in the continuum, it is hard to see any noncommutativity cleanly. Recently, Diirr 
and Hoelbling have looked instead at ml^ipip) / g"^ {g is the coupling) for staggered fermions 



22| | . They find clear evidence that the chiral and continuum limits do not commute for this 



quantity. In addition, the behavior is is very similar, qualitatively, to that observed for 

in the one-flavor Schwinger model. At the least, this shows that noncommutativity of limits 

is not inextricably tied to the rooting procedure in these 2D models. 

Finally, one may wonder whether there are order-of-limits problems in 4D that are directly 
analogous to those in the one-flavor 2D Schwinger model, and not of the kind treated above 
in SXPT. For example, in QCD with a single flavor, the value of (ipip) in the chiral limit is 
probably crucially dependent on having exact zero modes j^l, just as it is in 2D, leading 
to noncommuting chiral and continuum limits. In normal, multiflavor QCD, however, I 
know of no standard physical quantities that have this kind of sensitivity to the existence 
of exact zero modes. For example, the topological susceptibility seems to be well-behaved, 
both in SXPT j^J and in simulations 3]. This is not to say that the staggered discretization 
errors are negligible — indeed they are large. But there is no commutativity problem; 
the continuum suppression of the susceptibility seems to be present whether one takes the 
chiral or continuum limit first. Similarly, I do not expect any problem in the QCD rj' mass. 
Because of the anomaly, the rj' mass is generated through disconnected meson diagrams j^l • 
In the continuum limit, staggered fermions have the correct anomaly (for infinitesimal chiral 
rotations) and this should be adjusted correctly by the fourth-root procedure, since 

the anomaly is essentially perturbative. At finite lattice spacing, I expect that SXPT will 
correctly capture the modifications in the t]' physics due to staggered fermions, at least at 
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vacuum angle 9 = 0,^ just as ordinary XPT correctly captures t]' physics in the continuum 
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